a. 



LIBRARY 

RE'EARC; ■'r OOTS Dl\ 
N’V: ■ !<:• 

i,;_ . \ 



NPS5^79^qO^ 

^NAVA[ ‘P08TGRADUATE SCHOOL 

Monterey, California 




HEAVY TRAFFIC RESULTS 
FOR SINGLE SERVER QUEUES WITH 
DEPENDENT (EARMA) SERVICE 
AND INTERARRIVAL TIMES 

by 

P. A. Jacobs 
April 1979 



Approved for public release; distribution unlimited. 



FEDDOCS 
D 208.14/2: 
NPS-55-79-009 



NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 



Rear Admiral T, F, Dedman . J. R. Borsting 

Superintendent Provost 



Reproduction of all or part of this report is authorized. 
This report was prepared by: 



UNCLAS SIFIE D 



security classification of this page (When Data Entered) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


t REPORT NUMBER 

NPS55-79-009 


2. GOVT ACCESSION NO. 

oiC 


3. RECIPIENT'S CATALOG NUMBER 
7 


4 TITLE rend Subl/f/e; 

Heavy Traffic Results for Single Server Queues 
with Dependent (EARMA) Service and Interarrival 
Times 


S. TYPE OF REPORT a PERIOD COVERED 

Techn ical 


6. PERFORMING ORG. REPORT NUMBER 


7. AUTHORrs; 

P. A. Jacobs 


8 CONTRACT OR GRANT NUMBERrsj 


9 performing organization name and address 

Naval Postgraduate School 
Mon te re y , CA 9 39 40 


10 program Element project task^ 
AREA a WORK UNIT NUMBERS 


11 controlling office name and address 
Naval Postgraduate School 
Monterey, CA 93940 


12 report date 
April 1979 


13. number of pages 

27 


14 monitoring agency name a ADDRESSr^/ different from ControlUnfi Office) 


15 security class, (ot this report) 

Unclassified 


i5« declassification downgrading 
schedule 



16 Distribution statement (ot this Report) 



Approved for public release; distribution unlimited. 



17 distribution statement (ot the abstract entered in Biock 20. it ditterent from Report) 



18. SUPPL EMEnT ARY notes 



19 KEY WORDS (Continue on reverse aide if necessary and identify by block number) 

Cross correlated interarrival and service times 
Multivariate exponential random variables 
Heavy traffic 
Weak convergence 

EARMA sequences 

20. ABSTRACT (Continue on reverse side if necessary and identify by biock number) 

Models are given for sequences of correlated exponential interarrival and 
service times for a single server queue. These multivariate exponential 
models are formed as probabilistic linear combinations of sequences of inde- 
pendent exponential random variables and are easy to generate on a computer. 
Limiting results for customer waiting time under heavy traffic conditions 
are obtained for these queues. Heavy traffic results are useful for analyz- 
ing the effect of correlated interarrival and service times in queues on 



1473 EDITION OF 1 NOV 65 IS OBSOLETE 
S/N 0102-014-6 60 1 



1 JAN 73 



SECURITY CLASSIFICATION OF "^hiS PAGE 



Data Rntare>n 



unclassified 

_^uN|TY classification of this PAGETHhen Data Entered) 



20. ABSTRACT 

such quantities as queue length and customer waiting time. They can also 
be used to check simulation results. 



security classification of this PAGEfVVhon Data Entered ) 



HEAVY TRAFFIC RESULTS FOR SINGLE SERVER QUEUES 



WITH DEPENDENT (EARMA) SERVICE 
AND INTERARRIVAL TIMES 

by 

* 

P. A. Jacobs 

Naval Postgraduate School 
Monterey, California 

0 . ABSTRACT 

Models are given for sequences of correlated exponential 
interarrival and service times for a single server queue. These 
multivariate exponential models are formed as probabilistic linear 
combinations of sequences of independent exponential random vari- 
ables and are easy to generate on a computer. Limiting results 
for customer waiting time under heavy traffic conditions are 
obtained for these queues. Heavy traffic results are useful for 
analyzing the effect of correlated interarrival and service times 
in queues on such quantities as queue length and customer waiting 
time. They can also be used to check simulation results. 



Support from NSF grants ENG-77-09020 and ENG-79-01438 and 
ONR contract NR-42-284 is gratefully acknowledged. 



INTRODUCTION 
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Much of the work in queueing theory deals with models in 
which assumptions of independent service and interarrival times 
are made. One way to examine the effect of correlated interarrival 
and service times on measures of performance in queues is to 
consider models in which the interarrival and service times are 
correlated stationary sequences of random variables having given 
marginal distributions and compare results for these models with 
those for which the interarrival and service times are independent 
and identically distributed with the same marginal distributions. 

In [3] Jacobs and Lewis introduced a scheme for generating 

sequences of dependent exponential random variables. The 

EARMA (1,1) (exponential mixed autoregressive moving average with 

both autoregression and moving average of order 1) model is defined 

as follows. Let be a sequence of independent exponential 

random variables with positive finite mean v Let {J } and 

{K^} be independent sequences of independent {0,l}-random variables 

with P{j =1} =1-6 and P{K =1} = 1-p where 0 < 3 1 

n n — — 

and 0 < p < 1 are given constants. We define 



( 1 . 1 ) 

where 



X = 6 E + J A , , 

n n n n-1 



n = 1,2, 



( 1 . 2 ) 



A = pA T + K E 

^ r\ — \ -n in ' 



n = 1,2, 



Vie assume that Aq is an independent exponential random variable 
with mean v Under this assumption ^ ^ is a stationary 

sequence of dependent exponential random variables with mean \> ^ 
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The dependence is not Markovian in general and as was shown in [3] 

(1.3) Corr{X^,X^^^) = ( lp3 ) [ 6 ( !-(■ ) + (l-3)p] . 

Note that if 3 = 1 or 3 = 0 and p = 0, then is a 

sequence of independent exponential random variables. 

The EARMA(1,1) process is formed as a probabilistic linear 
combination of independent random variables. This technique can 
be extended to form sequences of multivariate exponential random 
variables. This idea was briefly explored in Lewis and Shedler [7]. 
These models have the advantage that the marginal distribution and 
correlation structure of the sequence are specified separately. 

In section 2 we will give several ways in which to generate 
sequences of dependent multivariate e.xponential random variables 
from probabilistic linear combinations of independent exponentials. 
These stationary multivariate EARMA sequences of dependent exponentials 
can be used as interarrival and service times in queues. Results 
for measures of performance for these queues can then be compared 
to queues with stationary independent sequences of interarrival 
and service times having the same marginal (exponential) distri- 
butions. Since the correlation structure in the multivariate 
EARMA sequences is specified by parameters independent of the 
marginal distribution, any difference between results for queues 
with these multivariate interarrival and service sequences and 
the corresponding queues with independent interarrival and service 
sequences having the same marginal distributions will be due to 
the correlation structure. 

Suppose that the interarrival and service times of a single 
s. rver queue form stationary multivariate EARMA sequences 
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of dependent exponential random variables. Results of Loynes [9] 
give conditions for the existence of limiting distributions for 
the waiting time of a customer and the queue length. However, it 
seems to be difficult to obtain exact results even in the simplest 
queues with EARMA interarrival and service times (cf. Jacobs [4]) . 

As a result we turn our attention to the possibilities 
of obtaining limiting results of the heavy traffic type for 
queues with stationary multivariate EARMA sequences of dependent 
interarrival and service times. Results of this kind will indi- 
cate how dependence in queues introduced by the EARMA scheme 
affects quantities such as customer waiting time and queue length, 
at least under heavy traffic conditions. These limiting results 
can also be used to check simulations. 

In this paper we will show that Kingman's [5] result 
for the equilibrium waiting time W in heavy traffic holds for 
queues having multivariate EARMA interarrival and service times; 
that is, we will show under conditions detailed in Section 3, that 
if the traffic intensity Ay ^ of the queue is close to 1, then 
the distribution of (1 - Ay ^)W is approximately exponential. 

The mean of this exponential distribution will in general be 
different from the case of independent interarrival and service 
times. Positive cross-correlation between the interarrival and 
service times will tend to decrease the mean; while positive 
correlation within the interarrival and service time sequences 
will tend to increase it. 
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In Section 3 we will state the heavy traffic result more 
precisely and present several examples of its application. Its 
formulation for multivariate EARMA queues is very dependent on 
the fact that the marginal distribution and correlation structure 
are specified independently of one another in EARMA formulations. 
In Section 4 we will give the proof of the heavy traffic result. 
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2. QUEUES WITH EARMA INTERARRIVAL AND SERVICE TIMES 



We will consider a single server queue at which customers 
arrive and are served in the order of their arrival. The 0th 
customer arrives at t = 0 and finds the server free. The 
arrival time, service time, and waiting time (excluding service) 
of the nth customer are denoted respectively by ~ ^k-1 ' 

S^, and W^ . In this section we will present some examples of 
the use of the EARMA scheme to model dependent interarrival and 
service times in queues. 

(2.1) EXAMPLE. Let {X^} be an EARMA (1,1) process with 
positive finite mean A ^ and parameters 0 < ^ 1 and 

0 < < 1; let EARMA (1,1) process independent of 

{X^} with positive finite mean p ^ and parameters ^ i. ^2 — ^ 
and ^ ^ P 2 ^ Note that if ^ ^ = 1 , = 1 , on = 1/ 

6 ^ = 0 / ^2 ~ = 0 ^ ^1 “ ^2 ~ = 0 , = 0 

B 2 = 0 » P 2 “ 0 ' then the queue reduces to an M/M/1 queue. 

Other special cases include queues with Poisson arrivals and 
EARMA (1,1) service times independent of the arrival process, 
and queues with EARMA (1,1) interarrival times and independent 
exponential service times. 

The sequence {(X^,S^)1 is stationary and 
{S^} are independent. However the ® (respectively S^^'s) 

are positively correlated with 
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(2.2) Cov(X^,X^^^) 



= A ^(l-6j^) (Sj^(l-Pj^) + (I-Sj^jpj^l 



(respectively 



(2.3) Cov(S^,S^^j^) = M ^P 2 ^ ( 1 - 62 ) [^2 (I-P 2 ) + (1-^2^ ^2^ ^ 



f rom (1.2). 



1 2 

Let (respectively F^) be the o-alqebra generated 

by Xj^,...,X^ (respectively Sj^,...,S^) and (respectively 

2 

G^^]^)' be the a-algebra generated by ' ^m+K+l'"'* (respectively 

®m+lA' Vk+r ■•■’• (respectively ) 

denote the collection of real-valued functions having finite 

second moment that are measurable with respect to F^ (respective!' 



=m 



G^,,), i = 1,2. If f. t L^(F^) and g. ^ L'^(G^ , ), i = 1,2, 
=m+k 1 =m =m+k 

then there is a k^ independent of f^, g^, i = 1, 2, such that 



2 ,^i 



‘0 



for k 2. 



( 2 . 4 ) |E[ (f^f 2 ) (g^L'^ 2 ^ ^ ■ ^f^i^ 2 ^ t^[qig 2 ]l 

= lElf^g^] E[f2g2] - E[f^] E(f2l E[g^] E[g2]l 

= |E[f^g^] (E[f2g2] - E[f2l E[g2]} 

+ E[f2] E[g2l {E[f^g^] - E[f^] E(g^]}| 
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(2.5) 



+ E[|f2l] E[|g2l] 5p^/^ E[f^]^/^ E[gJ]^/^ 
< 5[p^/^ + p^/^] E-[ (f^f^) E[ (g^g^) 



where (2.4) follows from the independent of ^^n^ 

and (2.5) follows from (6.5) of Jacobs and Lewis [3]. Hence, 

if f (respectively g) has finite second moment and is measur- 

1 2 

able with respect to the product o-algebra x (respectively 



=m 



1 '> 

G , X G , , ] 
=m+k =m+k 



then for k sufficiently large 



:2.6) |E[fg] - E[f] E[g]| £ 5E[f^]^/^ E[g^]^/^ 



Hence, the sequence {(X^,S^)} is asymptotically uncorrelated 
in the sense of Rosenblatt [10] and is <t)-mixing in the sense of 
Billingsley [1]. 



(2.7) EXAMPLE. In Lewis and Shedler [7] one of the schemes 

that is proposed for cross-correlating the arrival process and 

the service times is the following. Assume ^^n^ 

are independent sequences of independent exponential random 

variables with positive finite means A ^ and p ^ respectively. 

The nth interarrival time is assumed to be X and hence the 

n 

arrival process is Poisson with rate ■. Let the nth service time 

(2.8) S = HE + J (Au”^B ) 

n n n n 
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where 



(2.9) 



B = pB . + K X 

n n-1 n n 



and 3, Pt d-D and (1.2) so that 

(Sj^l and {^j^} cross-correlated sequences. We will 

assume that B^ has an exponential distribution with mean A 

Since E and B are independent, {S } is a stationary 

n n n 



sequence of exponential random variables with mean y 



-1 



If 



3 =1, then the queue reduces to an M/M/1 queue. The 
can be shown to be positively correlated with 



( 2 . 10 ) 



Cov(Sj^, ^n+k^ ^ ^ ^(1-3)^ k = 1 , 2 ,... . 



Further, the nth service time is positively correlated with 

the interarrival times X , X X, with 

n n-i 1 



( 2 . 11 ) 



Cov(S , X ,) = (Ay) ^(l-3)(l-p)p^ 
n n “K 



for k = 0 , . . . , n-1 . 

Note that (X,S,...,X,,,S,,) is a function of 
n n n+k n+k 

(B ,,X ,E ,J ,K ,X ,wE ,.,K X -,E ,J ^,,K .) 

n-1 n n n n n+1 n+1 n+1 n+1 n+k n+k n+k n+k 

Hence, since the sequence is stationary so is the 

sequence { 1. Further, by arguments similar to those 

used in obtaining (6.5) of Jacobs and Lewis [ 3], if f and 
g are as in (2.6) , then for k sufficiently large 



( 2 . 12 ) 



E[fg] - E[f] E[g]| < E[f^]^/^ E[g^]^/^ 



so that the sequence I (X^,S^) } is 4)-mixing. 
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There are many ways to cross-correlate the interarrival 
and service times using the EARMA scheme beyond those given 
above and in Lewis and Shedler [7]. We will present one more, 
somewhat more general than that of Example (2.7) in which the 
service times were coupled back over the previous interarrival 
times . 



(2.13) EXTUyiPLE. Let and {C^^} be independent 

sequences of independent exponential random variables having 
respective positive finite means A p ^ , and 1. Let 

X = 3E + J (A"^A ) 
n n n n 

S = yE + I (p~^A ) 
n ' n n ^ n 



(2.14) 

(2.15) 



with 



(2.16) 

whorG Pf 3 .nd S-TG as in (1.1) and (1.2) and 

is a sGquGncG of indGpGndGnt binary random variablGS indGpGndont 
of GvGrything with P{l^=l} = 1 - P{l^=0} = (l-y) for somG 
fixGd constant 0 1 Y <_ 1 • Again wg assumG that Aq is an 
indGpGndGnt random variablG having an GxponGntial distribution 
with mGan 1. Then {(E^/ ^n^ ^ ^ stationary 

SGquGncG of random variables and; hence, the same is true of 



{(Xn,S„)). 



9 





If r=l, B=1 or Y=l/ 3=0, p=0 or y = 0 , 


B = 1, 


p = 0, then the queue reduces to an M/M/1 queue. If 


B = 1, 


p = 0 then the queue reduces to a special case of Example 


(2.7) . 


Straightforward manipulation shows that in general the 



interarrival times are positively correlated with 



(2.17) 


Cov(X^, X^^j^) = A ^(l-B)^p^ , k = 1,2 



The service times are also positively correlated with 



(2.18) 


Cov(S^, = y ^(l-Y)^p’^, k = 1,2 



The interarrival times and service times are positively cross- 
correlated with 



(2.19) 


Cov(X^, = (1-3) (1 -y) (Ay) ^ p’^ , k = 0,l,... ; 



and 



(2 . 20) 


Cov(X^^j^, S^) - (1-y) (1 -y) (Ay) ^ p*^ , k = 0,l,... . 




Finally, by the arguments used in obtaining (6.5) in 



Jacobs and Lewis [3] , if f and q are as in (2.6) , then 
for k sufficiently larqe 

(2.21) |E[fg] - E[f] E[g]| ^ E[f^]^/^ E[g^]^/^ . 

so that the sequence ({)-mixing. 
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3. A HEAVY TRAFFIC RESULT FOR EARMA QUEUES 



Let (respectively ^he sequence of 

interarrival times (respectively service times) for a single 

server queue with first-in-first-out service discipline. 

Put U = S - X and assume that {U } is a stationary 
n n n n ^ 

sequence with E[U^] < 0. Let be the waiting time of the 

nth customer and W = sup . a )i i U, (with the convention 

^n 0 ^k=l k 

V 0 

that 2.^-1 “ ^ ' Loynes [9] has shown that if the sequence 

{U } is metrically transitive, then the distribution of W 
tends in the limit to the distribution function of W. Since 
it is usually difficult to obtain the distribution of W 
analytically, particularly when the service and interarrival 
times are dependent, we will obtain a heavy traffic limiting 
result for it. In this section we will give the statement 
of the result as well as some examples of its use leaving the 

proof to the next section. 

For each a > 0 let be a single server queue 

with FIFO queueing discipline having stationary sequences of 

interarrival times {X^ (a) } with positive finite mean A (a) 

and service times {S (a)} with positive finite mean u (ut) ^ 

n 

such that, if U (a) = S (a) - X (a), then {U (li)} is a 

' n n n n 

stationary sequence of random variables. Let 

(3.1) E[U (a)] = v(a) < 0; Var[U (a)] = s^(a) ; 

n n 
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(3.2) 



(3.2) 

(3.3) 


oo 

a^(a) = Var U, (a) +2 1 Cov (U^ (a) (a) ) ; 

n=l 

_2 

d(a) = |v(ot)|/a(a) and n(a) = d(a) 



We will assume the following: 

(3.4) ASSUMPTION. There exists a distribution function G so 

2 

that lim P{[U (a) - v(a)] < x} = G(x) for all continuity 

a ^ n — 



points 


2 2 

X of G; lim v (a) = 0 ; and 0 < lim s ( x) = s < . 

a 


(3.5) 


ASSUMPTION. 0 < lim ^^( J = < °° . 

j -> • 


(3.6) 


ASSUMPTION. If f (respectively g) is a function with 


finite 


second moment measurable with respect to the a-algebra 



generated by ( a) , . . . , ( a) (respectively U^+]^(a), ( a) , . . . ) , 



then 




(3.7) 


|E[fg] - E[f] E[g]| ^ c{)(k) E[f^]^/^ E[g^]^/^ 



for some decreasing function (p bounded above by 1 not depending 



on a 


satisfying 


(3.8) 


00 

1 4>(n) < « . 

n=l 



Assumption (3.4) concerns the convergence of the distribution of (a) 
_ 2 

- ^ first two moments as a Assumption 

(3.6) is roughly a condition that the sequences {U^(a), n=l,2,...} 
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be 4)-mixing uniformly over all a. This assumption guarantees 
that the series in (3.2) converges and 



lim 

n oo 




Assumption (3.5) concerns the convergence of this normalized 
variance term as a . 



where [x] is the largest integer less than or equal to x. 

Let ID(I0,“)) be the collection of real-valued, right continuous 
functions on [O,"’) which have left limits everywhere endowed 
with the metric given in Lindvall [8]. 

(3.9) THEOREM. If (3.4)-(3.6) are satisfied, then the process 
= {Z^(t);t 2l 0} converges weakly in JD([0 ) ) to a 

Brownian motion with negative unit drift. 

The proof of this result will be given in the next section. 
Le t 



Put 



[n (a) t ] 



Z (t) = d (a) o (a) 
a 



t > 0 



n 

W(a) = sup ^ 
n > 0 k = l 



I • 



Then it follows from Theorem (3.9) that 
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a 



a 



lim 

-> CO 



p 



fd(a) 

lo(^0 



W(a) 



< 




lim P{ sup Z (t) x} 
-> >» t > 0 ^ 



= 1 - exp{-2x} 

where the last equality follows from the continuous mapping 
theorem (of. Billingsley [1]). Thus we have the following 

resul t . 



(3.10) THEOREM. If (3.4)-(3.6) hold, then 

"" 2 

lim P{ U ( j) 1 W(a) < x} = 1 - exp{-2o x) • 

X ^ " 



Hence, if the traffic 

than but close to 1 and 

ing (3.7) and (3.8), then the 

1 2 

exponential with mean “ '2 ^ 



intensity Xy of a queue is less 
is a stationary sequence satisfy- 
distribution of W is approximately 
E[Uj^] ^ where 



= Var 






-f 2 



I 

n=l 



Cov(U, ,U 



1+n' 



Thus positive cross correlation between the interarrival and service 
times will tend to decrease the average mean waiting time- while 
positive correlation within the interarrival time or service 
time sequence will tend to increase the mean waiting time. 
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(2.11) EXAMPLE. For each a > 0, let be a queue as in 

Example (2.1) with {X^(a)l (respectively {S^(a)}) forming 
an EARMA(1,1) process with parameters X (a) , 3^^ and 
(respectively pi (a), and P2) • Assume A (a) < pi (a) for 

each a, but lim X(a) = lim , (..) = m; then hypotheses 

a 06-^ 

(3.4)-(3.6) are satisfied and Theorem (3.10) holds with 

= 2 m"^{l + ( 1 -pj^)"^ ( 1 - 63 ^) [ (l-Pj^) 63 ^ + ( 1 - 33 ^)P 3 _] 

+ (1-P2)"^ ( 1 - 32 ) [ (I-P 2 ) ^2 + (l-32)p2n . 

Hence, if the traffic intensity Api * of a queue with independent 

EARMA(1,1) services and interarrival times is less than but 

close to 1 , then the distribution of the equilibrium waiting 

1 2 

time W is approximately exponential with mean ~ 2 ^ £[ 113 ^] 

-1 "1 2 

where E[U 3 ^] = p - A and o is given by 

0^ = A"^ {1 + 2(1-P3^)"^ (1-33^) [(1-P3_)33^ + (l-33^)p2]) 

+ p"^{l + 2(1 -p 2)"^ (l-e2^p2^^ • 

Thus the dependence in this queue tends to increase the mean 
customer waiting time at least in close to heavy traffic 



conditions . 



(3.12) EXAMPLE. For each a > 0^ let Q be a queue as in 
Example (2.7) with parameters E[X^(a)] = A(a)“^, 3, and p 
and E[S^(a)] = y (a) Again we will assume A (a) < y(a) 

for each a but lim A(a) = A = lim y(a) = y = m. Hence 

(Jt“ ^ oo 

(3.4)-(3.6) are satisfied and Theorem (3.10) holds with 

= a"^ + y“^ - 2(1-3) (1 - P) (Ay)“^ 

+ 2p(l-p)"^ [(l-3)^y“^ - (1-3) (1 - p) (Ay) ~1] 

= A ^ + y“^[l + 2(1-3)^ p(i_p)"l] _2(Ay)“^ (1-3) 

= 2m"^[l + (1-3)^ p(l-p)"^ - (1-3)] . 

2 2 

The quantity a - can be either positive or negative where 

^ = A ^ + ). ^ = 2m ^ . It is negative if p < (2-3) ^ 

ind 

and nonnegative otherwise. Hence, if p is sufficiently less 
than (2-3) then the average customer waiting time in the 

queue will tend to be less than the waiting time in the inde- 
pendent case at least in close to heavy traffic conditions. 

It will tend to be greater if p is sufficiently greater than 

(2-3) 

(3.13) EXAMPLE. For each t > 0, let be a queue as in 

Example (2.13) with parameters A(u), y(a), p, 3 , and y. 

Assume lim A( .;) = A = lim y(a) = y = m as before. 

O '- ' a->- °° 

Then (3.4)-(3.6) are satisfied and Theorem (3.10) holds with 
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- 2(Xp) 



(1-6) (1-Y) 



2 

a 



-2 , ^ -2 

]j + A 



-1 



+ 2p(l-p) ^(1 -y) -a ^(1-6)}^ 



= 2m ^ [l-(l-6) (1 -y) + P(l-P) ^{(6 -y)}^] . 



2 2 

The difference a - can be either positive or negative; 

2 -2 -2 -2 

(a. , = A + u = 2m ) . The difference will be negative if 

ind ^ 

2 

p < (1-6) (1-y)/[(1“3) (1“Y) + (6“y) ] nonnegative otherwise. 

Hence under close to heavy traffic conditions, the average 
customer waiting time will tend to be less than the independent 
case if (1-6) (1-y)/[(1~B) (1-y) + (6-y)^] is sufficiently greater 

than p. It will be greater than the independent case if the 
same quantity is sufficiently less than p. 
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PROOF OF THE MAIN RESULT 



4 . 

In this section we will give a proof of Theorem (3.9) . 
We will assume that (3.4) -(3.6) are satisfied. Let 



Y.. (t) 



d (a) 
a (a) 



[n (a) t] 

I 



k=l 



[\ (a) 



v(a)] . 



By Theorem 3' of Lindvall [8] and the continuous mapping theorem. 
Theorem (3.9) will be proved if for each a > 0 we can show 
the following result . 

(4.1) THEOREM. The process {Y^(t); 0 < t < a} converges 
weakly on ID([0,a]) to Brownian motion. 

0D([O,aj) is the space of real-valued right continuous functions 
on [0,a] with left-hand limits endowed with the Skorohod 
topology . ) 

We will prove Theorem (4.1) for a = 1 by a series of 
lemmas. The proof for arbitrary a > 0 is similar. 

(4.2) LEMMA. The process {Y^(t) ; 0 < t < 1} has asymptotically 
independent increments. 

Proof. Let 0<s. <t. <ST<t„<'""<s <t <1. Let 

-1—1 2—2 r— r — 

C. be the event {Y (t.) - Y (s-) < x.} for real x.. Then 

1 <1 a 1 — 1 1 

is an event in the o-algebra generated by 
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U[n(a)s.]+1' ••• ' ^[n(a)t.]- smallest difference 

- t^_^, then [n(cx)s^] + 1 - [n(a)t^_j^] >_ [n(a)6]. By (3.7) 

r r 

|P( II c. ) - n P(C. ) I < r(|)( [n(a) 6] ) . 
i=l ^ i=l ^ “ 

Since 6 > 0, rc|) { [n (a ) o ] ) tends to zero as a ^ and the 
result follows. 

In the following results, let D^(a) = U^(a) - v (a) 

r-m 

and V (a) = V. , D. (a) . By (3.4)-(3.6) the following result 
m ^1=1 1 

holds . 



(4.3) LEMMA. 



CO 



a) ^E[Vj^(a)^] ^2 |E[D^(u) D^^.(,;x)]|; 

b) lim n(a) ^ E[V^, . (a)^] = ; 

Qt -► OO 



C) 



If is bounded by 

where depends on 

<t> 



C, then E[V , . 

n (a) 

Ip alone. 



4 4 

(a) ] < K C 

- P 



n (a) 



2 



(4.4) LEMMA. The sequence of random variables 
-1 2 

{n (a) (a) ^ ; a > O) is uniformly integrable. 

2 

Proof . By (3.4), {Dj^(a);a^0l is uniformly integrable. The 

remainder of the proof is very similar to that of (20.48) on 
page 176 of Billingsley [1] and will be omitted. 
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(4.5) LEMMA. For each positive e there exists a y, y > o 



and an integer 



n^ such that n(a) n^ implies 

P{ max |v. (a) I >_ y/n (a) } < . 

i < n (a) ^ y 



^ * 

Proof . Put Cj^(a) = |Dj(a)|. By (3.4), there exists 

an increasing sequence of integers m^ such that n (a) m^ 

2 2 —1 

implies P{|Dj^(a)| y /n (a) /i } £ (y n(a)i ) for each 
positive f . Let p(u) = i if n(a) < + 2 ' then p(a) 

tends to infinity as a. ^ but 

(4.6) lim n(a) P{C , .(a) >y /n (a) } 

a -> x> 

_< lim n(a) p(a) P{|Oj^(a) | >_ y /n (a) p(a)~^} 

a 

_< lim n(a) p(a) P{|D^(a) | >_ y /n (a) p(a)~^} 

a -► oo 

2 -1 

^ lim n(a) p (rt) [y n (a) p(a) ] = 0 . 

■I -y ■ -r; 

We can further choose p(a) so that p(a) ^ n(a). Fix e > 0. 
By the previous lemma we can choose y so that y > a and 
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(4.7) ^ ^ < £Y 

for all a. Let 



then 



A . ( a ) = { max I V . (a ) I < 3y /n (a) < I V . ( a ) | } ; 
^ j < i ^ 



P{ max |v. (a) I >_ 3y /n ( a ) } 
i < n (a) ^ 



<_ P{ I Vr (a) I >_ Y /n (a) } 
n (a) -1 

+ J P(A^(a) II - V^(a)| > 2y /nToT} ) . 



The summation is bounded above by 



n(a)-p(a)-l 

I P{|v^(a) - '^i+p(„) I /nTaD 



n (a) -p (a) -1 

+ P(A.(a) n Uv_,(^,(a) - V.^p,^,(a)| >Y /iUCT) ) 

n (a) -1 

+ I P{|V , . (a) - V. (a) I >_ Y /n (a) } 

i=n (a) -p (a) ^ ^ 



Each term in the first and third of these sums is bounded above 
* 

by 2 . 7 /n (a) } . Since A^ (a) is measurable with 

respect to the a-algebra generated by U^(v) by (3.7) 
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p{ max |v. (a) I >_ 3y /n ( a ) } 
i ^ n (a) ^ 

_< P{ IVr (^) (a) I >_ Y >/n (a) } 

+ [n(a) + p(a) ] (“) 2L (a) } 

n (a) -p (a) -1 

+ \ P(A^ (a) [P{ (a) - (a) I > Y /nT^}+(J) (p ( u) ] 

^ 2£Y~^+ct' (p (-. ) ) + [n(a) +p(a)] P{Cp^^j(a) >_ y /n ( a ) } 

by (4.7) and the fact that the A^(a) are disjoint. Hence by 
(4.6) and (3.8) for a large 

P{ max |v. (a) | 3y /n (a) } ^ 3 ey~^ . 
i ^ n (a) ^ 

PROOF OF THEOREM (4.1) . For a = 1: By Lemma (4.5) the 

sequence tight on ]D[0,1]. Further lim E[Y^(t)] = 0 

and 

( a) 

by (3.5). Since Y,^(t)^ = [^^(0 n(a)]"^ ^ [n ( a) t ] ^ ^ ' 

{Y^(t)} is uniformly integrable for each t by Lemma (4.4). 

By Lemma (4.2) {Y^(t) ; 0 < t < 1} has asymptotically independent 

increments. The result now follows from Theorem (19.2) of 
Billingsley (!]• 




lim E[Y (t)^] = lim ( /n ( a ) o(a)) ^ E 

^ a 

^->■00 (^->*00 



[n in) t ] 



/ in wx; 
' k = l 



D 
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EXTENSIONS 



5 . 

This paper presented a heavy traffic result for the 
limiting distribution of the customer waiting time in single 
server queues with dependent interarrival and service times. 
There are many different ways in which to model dependent 
interarrival and service times in queues using the EARMA scheme. 
Some other ways are given in Lewis and Shedler [7]. Properties 
of these multivariate exponential models will be given 
elsewhere . 

Work is currently going on in the simulation of these 
queues. Some simulation results have been reported for the 
queue of Example (2.7) with B = 0.5, p = 4 and various 
values of P and A in Boonsong [2] . They seem, to indicate 
that the heavy traffic result approximates the average 
customer waiting time fairly well for a traffic intensity of 
.9 and quite well for a traffic intensity of .95. However, 
the distribution of W in a simulation with p = 0.25 and 
traffic intensity .995 shows a large amount of underdispersion 
relative to an exponential distribution. Thus, it seems that 
the convergence of the waiting time distribution is very slow. 

It is expected that further simulations will give 
more indication of the validity of using the heavy traffic 
result to approximate the average customer waiting time, or 
quantiles of the waiting time distribution under less than 
heavy traffic conditions. These results will be reported 
on elsewhere. 
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